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^ ■ Abstract 

Consider the Schrodinger operator —V + q with a smooth compactly 
supported potential q, q = q{x),x G R'^. 

Let A{I3, a, k) be the corresponding scattering amphtude, k'^ be the 
energy, a G S*^ be the incident direction, /? G 5^ be the direction of 
scattered wave, be the unit sphere in R^. Assume that k = ko > 
is fixed, and a = oq is fixed. Then the scattering data are A(/3) = 
^ . A{f3, ao, A;o) = Aq{f3) is a function on S'^. The foUowing inverse scatter- 

O I ing problem is studied: 

IP: Given an arbitrary f G L?'{S'^) and an arbitrary small number e > 0, 
can one find q G C^{D), where D G R^ is an arbitrary fixed domain, 
O ■ that \\Aq{P) - /(/3)||i2(52) < e? 

\ A positive answer to this question is given. A method for constructing 

such a g is proposed. There are infinitely many such g, not necessarily 
real-valued. 

H . MSG: 35R30; 35J10; 81Q05 

■ 
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1 Introduction 

Consider the scattering problem: 
Find the solution to the equation 

[V^ + A;^ - q{x)]u = in R^ (1) 

such that 



M = e*'="-^ + A(/3,a,fc)— + o( - ) , r = |a;|^oo, - = /3, (2) 
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where a G S*^ is a given unit vector, S*^ is the unit sphere, k = const > 0, fc^ 
is the energy, a is the direction of the incident plane wave uq := e*'^°''^ is the 
incident plane wave, /3 is the direction of the scattered wave. The function 

A{(3,a,k) = Ag{(3,a,k) (3) 

is called the scattering amplitude corresponding to the potential q{x). 

If q{x) G C^(R^) and is a real valued function, then the scattering problem 
([T])-(l2]) has a unique solution, the scattering solution. There is a large literature 
on this topic, see, for example, [T] and and references therein. The scattering 
theory has been developed for much larger classes of potentials, not necessarily 
smooth and compactly supported. 

We prove existence and uniqueness of the scattering solution assuming 
that Img < 0, see Lemma 2.5 in Section 2. The inverse scattering problem 
consists in finding q{x) in a certain class of potentials from the knowledge 
of the scattering data A{f3,a,k) on some subsets of the set S'^ x S'^ x R_|_, 
where R+ = [0, oo). If A{P, a, k) is known everywhere in the above set, then 
the inverse scattering problem is easily seen to be uniquely solvable in the 
class of Lq(R^), that is, in the class of compactly supported square-integrable 
potentials, and in much larger class of potentials. If the scattering data is 
given at a fixed energy, and A; = /cq > for all /3 G 5*^ and all a G S*^, then 
uniqueness of the solution to inverse scattering problem was proved originally 
in [2]. An algorithm for finding q{x) from the exact fixed-energy scattering 
data and from noisy fixed-energy scattering data was given in [3] , where the 
error estimates of the proposed solution were also obtained, see Chapter 5 in 

Only recently the uniqueness of the solution to inverse scattering problem 
with non- over- determined data. A{—(3, (3, k) and A{(3, ao, k) was proved, see [5], 
and [6], [9]. 

The data y4(— k),yp G S*^ and all > are the back-scattering data, 
the A{(3, aQ, k),\/f3 G and VA; > are the fixed incident direction data. The 
scattering data are called "non-over-determined" if these data depend on the 
same number of variables as the unknown potential, that is on three variables 
in the above problems. 

Note that the data A{(3, a, ko),Wa, (3 E S'^ and a fixed k = kQ > are over- 
determined: they depend on four variables while q depends on three variables. 

The inverse problem IP with the data A{I3) := Ag{f3) := /l(/3,ao;^o) is 
under- determined: its data depend on two variables. This problem, in general, 
does not have a unique solution in sharp contrast to the inverse scattering 
problems mentioned above. The IP was not studied in the literature. 

In this paper the IP is studied. Assume that D C R^ is an arbitrary fixed 
bounded domain. 

Let us formulate the inverse problem: 
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IP: Given an arbitrary f G L?'{S'^) and an arbitrary small number e > 0, 
find a q E C^{D) such that 

||A(/3)-/(/3)||^.(5.)<e. (4) 



The IP's formulation differs from the formulation of the inverse scattering 
problems discussed earlier: 

i) there is no unique solution for the problem we are discussing, 

ii) in place of the exact (or noisy) scattering data a function on 5*^ is 
given, which, in general, is not a scattering amplitude at a fixed a = Oq and a 
fixed k = ko corresponding to any potential from L^{D). 

The main results of this paper include: 

a) A proof of the existence of g G C^{D) such that (jlj) holds; 

b) A method for finding a potential q G C^{D) for which (jll) holds; 

c) an analytical formula for a function h = qu, where u is the scattering 
solution at a fixed k and a fixed a, corresponding to q. 

In section 2 we prove that the set {A[(])} corresponding to all q G C^{D), 
is dense in L^(5'^), and that the set of functions {h} = {qu} is dense in 
L?'{D) when q runs through all of L'^{D). Here u = u{x) = u{x, ao, fco) is the 
scattering solution corresponding to the potential q, that is, the solution to 
the scattering problem ([I])-© with a = ao E and k = k^ > 0. 

In section 3 an analytical formula for q is given. The q computed by this 
formula generates A{f3) satisfying (j4]). 

We do not discuss in this paper the relation of our results with the theory 
of creating materials with a desired refraction coefficient, see [Z], [H]- 

2 The density of the set in L^S'^) 

Let us start by proving Lemma the following lemma. 

Lemma 2.1 If the set {A{P)} is dense in L'^{S^) when q runs through all 
of L^{D), then it is dense in L'^{S^) when q runs through C^{D). 

Proof. The set C^{D) is dense in L'^{D) (for example, in L?{D) norm), and 
the scattering amplitude A{I3) depends continuously in the norm of L'^{S'^) on 
g, that is, 

(5) 

Estimate ^ follows, for example, from the known lemma of the author (see 
m, P-262): 

- 4:Tx[Aq^{(5) - Aq^{l3)] = I [qi{x) - q2{x)]ui{x,aQ,kQ)u2{x,-(3,kQ)dx, (6) 

JD 
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and the well-known estimate 

sup\u{x,a,k)\ < c, (7) 

where c > is a constant depending on the L^{D) norm of q and uniform with 
respect to A; G [a, oo), a > is a constant and a G S'^. Lemma 2.1. is proved. 
□ 

Thus, in what follows it is sufficient to establish the density of the set {Aq{l3)} 
in L'^{S'^) when q runs through L?'{D). 

Theorem 2.2 For any f G LF'{S'^) and any e > there exists a q & C^{D) 
such that estimate (jlj) holds, where A{P) = Aq{l3) is the scattering amplitude 
corresponding to q, and a = ao & S"^ , k = ko > are fixed. 

Proof. By Lemma [2.11 it is sufficient to prove that the set {y4(/3)} is dense 
in L^(S'^) when q runs through all of LF'{D). Assuming the contrary, one finds 
a function / G L'^{S'^) which is orthogonal in L'^{S'^) to any function A{(3). It 
is well-known that 

- AnAiP) = [ e~'^P yh{y)dy, h{y) := q{y)u{y), (8) 

J D 

where u{y) = u{y, ao, k^) is the scattering solution. Using the orthogonality of 
/ to A[(5), one gets 

I fW) [ e-"'^-yh{y)dyd(3 = 0, ^q e L\D). (9) 

We prove later that when q runs through all of L'^{D), the corresponding h 
runs through a set L dense in L'^{D). Taking this for granted, one can replace 
in ([9]) the expression Wq G L'^{D) by \/h G L'^{D), and rewrite ([9]) as follows: 

/ hiy) [ f{P)e-'^^-ydpdy = 0, V/i G L. (10) 

This implies 

/ f{P)e-'^P-yd(5 = Q, ^yeL'iD), (11) 

where k = ko is fixed. In what follows, we write everywhere k for ko and a 
for ao- The integral in (fTTj) can be considered as the Fourier transform of a 
compactly supported distribution 

^7(A,/3):=^(0:=/(|)^^|^, (12) 
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where ^ G R^, A = |^|, /3 = f , and 
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~9{y) ■■= giOe-'^-'d^ = X'dX e-''^yg{X,^)d^. (13) 
Jr^ Jo Js^ 

Since distribution ( IT2l) is supported on the sphere |^| = k, which is a compact 
set in R^, its Fourier transform is an entire function of y. This function vanishes 
in an open in R^ set D by f fTT]) . Therefore, it vanishes everywhere in R^. By 
the injectivity of the Fourier transform one concludes that = 0. Therefore 
the assumption that the set {y4(/3)} is not dense in L^(S'^) is false. Theorem 
12.21 is proved under the assumption that the set {h} is dense in LP'{D) when q 
runs through all of L'^{D). In Theorem 12. 4[ see below, this density statement 
is proved. Thus, one can consider Theorem 12.21 proved. □ 

Remark 2.3 If one defines A{f3) = — ^ f^e~'^'''^'^h{y)dy and assumes that 
{h} runs through a dense subset of L'^{D), then the corresponding set {A{/3)} 
is dense in L?'{S'^), as follows from our proof of Theorem \2.'A 

Theorem 2.4 The set {h{x)} is dense in L'^{D) when q{x) runs through 
a dense subset of L^{D), where h{x) = q{x)u{x), and u{x) is the scattering 
solution. 

Proof. If u is the scattering solution, then 

u{x) = Uo{x)- g{x,y)h{y)dy, Uo{x) := e''"^"' , h = qu, (14) 
Jd 

^ik\x-y\ 

9{x,y) = —. 7- (15) 

4:7r\x — y\ 

u{x) uo{x) - g{x, y)h{y)dy 

If the function q{x), defined in (fT6|) . belongs to L'^{D), then the function u{x), 
defined in (fT^ . is the scattering solution, corresponding to g G L'^{D), defined 
by formula (fT6ll . Uniqueness of the scattering solution is guaranteed if Img < 
by the following Lemma [2.51 

Lemma 2.5 Assume that Imq < and q G L'^{D), q = in D := R^ \ D. 

Then there exists a unique scattering solution, that is, the solution to problem 
dl])-©. This solution is also the unique solution to equation (IT^ if h = qu. 

Proof. It is sufficient to prove uniqueness of the scattering solution. Indeed, 
the scattering solution solves a Fredholm-type Lippmann-Schwinger equation, 

u{x) = Uo{x)- / g{x,y)q{y)u{y)dy, 
Jd 



Define 
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and the uniqueness of the solution to this equation imphes the existence of 
this solution by the Fredholm alternative. 

Suppose that there are two scattering solutions, Ui and U2, that is, solution 
to problem ([I])-© (with a = and k = ko > 0). Then the function v := 
Ui — U2 solves equation ([1]) and satisfies the radiation condition at infinity: 

[V^ + P - qix)]v = in (17) 

— zkv = o - , r := X — 7- 00. (18) 

or \r J 

Multiply equation (fT7|) by v, the bar stands for complex conjugate, and the 
complex conjugate of (ITTl) by v and subtract from the first equation the second. 
The result is 

v{V^ - k^)v - v{V^ - k^)v - (g - q)\v\'^ = 0. (19) 

Integrate (HM over a ball Br of large radius R, centered at the origin, and use 
the Green's formula to get 

f I Ov dv \ f 

/ {v—-v—]ds-2i lmq(x)\v(x)\'^dx = 0. (20) 
J\x\=R \ dr drj Jb^ 

Using the radiation condition f[T5]) one rewrites (120]) as 



2ik 



I \v\^ds-2il \mq\v\^dx + o{l) = (21) 

J\x\=R Jbr 

where o(l) — )■ as i? — > 00. Thus, if Img < relation 0211) implies that 

lim / \v\'^ds = 0. (22) 



x\=R 



The radiation condition (fT8|) and equation (fT7|) with a compactly supported 
q G L'^{D), implies that f = in R'^, see Lemma [2. II on p. 25 in [TO]. 

The scattering solution solves equation f|T^ . because it satisfies the equa- 
tion 

u = Uo- / g{x,y)q{y)u{y)dy, (23) 
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SO that with the notation h := qu one gets from fl23|) equation ([T] 

Conversely, assume that u is defined by equation ([H]) and h in f[T^ is 
equal to qu, where q is defined by formula (1161) and q G L'^{D) . Then this u 
solves equation fl23l) . and, therefore, it satisfies equations ([1]) and ([2]). To check 
this, apply the operator + k"^ to equation (|T4l) and use the known formula 
(V^ + k'^)g = -6{x - y). 

The result is (V^ + k'^)u = h = qu, so equation ([T]) holds. The radiation 
condition holds because q is compactly supported and g satisfies the radiation 
condition. Lemma [2.51 is proved. □ 
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Although the scattering solution, in general, is not unique if Img > and may 
not exist in this case, it does exist, even if Img > 0, if g, defined by formula 
f lT6|) . belongs to L'^{D). Indeed, then qu = h and 

u{x) = uo{x) - / g{x,y)q{y)u{y)dy, 
Jd 

so that u := Uo{x) — g{x,y)h{y)dy is the scattering solution. 

Let us assume now that the function (fT6l) does not belong to L'^{D). Since 
h G L'^{D) the function u{x) G Hf^^{K^), where H^^^CR^) is the Sobolev space. 

Changing h slightly one may assume that his a smooth bounded function in 
D. Such a change leads to a small change of the scattering data A{(3) in L'^{S'^) 
norm. Thus, we will assume below that h is bounded in D in absolute value. 
The function f|T6l) does not belong to L'^{D) if and only if the denominator in 
ffTBl) has zeros. Let 

N := N{u) := {x : u{x) = 0,x e D}, Ns := {x : \u{x)\ <6,x e D}, 

where 5 > is a small number, and let Ds := D \ Ns. 

The idea of the argument below is to show that the set N is generically a 
line in R^, and that there exists a function 



such that 

\\h5-h\\L2(^D)<c5'. (25) 

Moreover, the function hs G L°°{D). This function can be made smooth by an 
approximation by a C^(D)-function. The corresponding smooth generates 
As{l3) which differs slightly from the original /(/?)• The corresponding ^^(x) 
can be defined as follows: 



U5{x) := uo{x) - / g{x,y)hs{y)dy, (26) 
Jd 

where 

'-={o^ '''' 

Consequently, we will prove that a small change of h may be arranged in such 
a way that the corresponding change of q leads to a potential qs which belongs 
to L°°{D). 



Lemma 2.6 The set N is a line in R^. 
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Proof. Let u = ui + iu2, where ui = Keu and U2 = '^u. Then the set is 
defined by two equations in R^: 

ui{x) = 0, ^2(0;) = 0, X e (28) 

The functions Uj e Hf^^(Il^), j = 1,2, because u G Hf^^(Ii^). Therefore, each 
of the two equations in ( 1251) is an equation of a surface. 

The two simultaneous equations fj28|) generically describe a line i in R^. 
By a small perturbation of h one may ensure that the line i := {x : Ui{x) = 
0,ui{x) = 0,x E D} is smooth in D and the vectors Vuj,j = 1,2, on i are 
linearly independent. 

Lemma 12.61 is proved. □ 

Lemma 2.7 There exists a function ([21]) such that (125!) holds. 

Proof. Consider a tubular neighborhood of the line i in D. This neighbor- 
hood is described by the inequality p{x, t) < 5, where p{x, t) is the distance 
between x and £. Choose the origin O on £ and lets the coordinates Sj,j = 1, 2, 
in the plane orthogonal to i, be directed along the vectors Vuj\i respectively, 
while the third coordinate S3 be directed along the tangent line to i. The 
Jacobian J of the transformation {xi,X2,X3) — )■ (si, 52,33) is non-singular, 
I J| -|- I J^^l < c, because Vuj,j = 1, 2 are linearly independent. 

Define hs by formula ( 124|) . us by formula ( 126|) . and qs by formula (127|) . 

Note that 

us{x) = u{x)+ g{x,y)[h{y) - hs{y)]dy. (29) 



JD 

Since hs = h in Ds and /i^ = in A''^, one obtains inequality 

Whs - h\\L2{D) = ||^||L2(iV,-) < c6^, 

because h is bounded and 5^ is proportional to the area of the cross section of 
the tubular neighborhood. Furthermore, 

\us{x)\ > \u{x)\ — c I — j -, c = max |/i(x)|. (30) 

Jjq^A'K\x-y\ xeNs 

Denote 

:= sup / (31) 
x&Ds Jns 4:7i\x-y\ 

By construction 

\u{x)\>5 if xeDs. (32) 
Therefore, inequality (l30l) implies 

\us{x)\>6 - I{6) yxeDs. (33) 
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Let us estimate as 5 — ?■ with the aim to prove that 

inf \us(x)\ > b(6) > 0, hm ^ = 1. (34) 

So, for sufficiently small 5 > one has b{5) > |. 

J holds, then |u5(x)| is strictly positive in Ds and, therefore, qs := ^ 
is a bounded function in Ds- 



If 

i b 

Since hs 

\qs - qllv^iDs) 



h in Ds, one obtains 
hs h 

Us 



u 



< \\h\ 



LHDs) 



U — Us 



UsU 



<c51n-. (35) 





Here the following inequalities for x & Ds were used 

5^ 



C{Ds) < C, 

usix) — uix)\ < 



usu\ > — in Ds, 



\g{x,y)h{y)\dy < c5 



dy 



4:Tt\x — y\ 



< cSI(S). (36) 



Here and throughout by c > various constants, independent of 6, are denoted. 
To estimate I{S) one argues as follows: 



dy 



HS) <— ri<c pdp 



4vr Jn, \y 



c5 



dS3 



'P' + si 



<c5Hn\, 





(37) 



where the unit 1 is a finite coordinate along the S3 axis and we have used the 
following estimate: 



ds'i 



'P' + si 



In ( S3 + \Ip^ + si 



' , 1 
< cm-, 

P 



(3^ 



where c > 1 is a constant. Estimates ( ]33l) and (1381) imply inequality (E 
Thus, the existence of a function fl24l) is proved, and 



\hs 



L'^Ns) 



< c I dx < cS . 

Ns 



(39) 
□ 



Lemma 12.71 is proved. 

From Lemma [2.71 the conclusion of Theorem 12.41 follows. Indeed, if /i G L'^{D) 
is an arbitrary function and q, defined by formula f|T6l) . belongs to L^{D), then, 
as was proved above, the u, defined by formula f|T^ . is the scattering solution. 
This scattering solution is unique if Img < by Lemma 12.51 If q, defined 
by formula flTBl) . does not belong to L'^{D), then there is a bounded function 
qs{x), approximating q, see fl5S]) . such that the corresponding hs approximates 
h{x) well, see (!39l) . and the corresponding us{x) is the scattering solution 
corresponding to qs{x). 

Theorem 12.41 is proved. □ 
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3 Formulas for solving inverse scattering prob- 
lem with fixed a and /c > 

The inverse problem (IP) was formulated in the Introduction. Given e > 
and an arbitrary G L'^{S'^) we first find h{x) G L'^{D) such that 



fi(3) + -i- / e-^'Pyh{y)dy 



D 



< e. 



(40) 



This can be done (non- uniquely !) in many ways. Let us describe one of the 
ways. Without loss of generalities assume that D = B = is a ball of radius 
R centered at the origin. Expand the plane wave e~*'^^'^ and h{y) into the 
spherical harmonics series: 



-ik(5-y 



J2^^^'Ukr)Yii-/3)Yi{yO) 



r = \y\, y 



._ y 



(41) 



1=0 



/ vr \ 1/2 
ji{r) := I — ) J,^i(r 



An 



\2rJ 



21 



- J2 YUx')YUy')=Pi{x'-y'), (42) 



m=—l 



Ji+i{r) is the Bessel function regular at the origin, Yi{a) are the spherical 



harmonics: 



Yi{a) = Yi^rnia) = 
— I < m < I, 
PirnicosB) = (sin( 



-1) 2 'i' 



/47r 



(2/ + !)(/- 


m\)\ 


{1 + 


m\)\ 



1/2 



1 d\t^-l 



2H\ dt^ 



e'^^P^cose), 

(43) 
(44) 



t = cos6, the unit vector a is described by the spherical coordinates {6, ip), < 
ip < 27T, 0<6<7T,—l<'m<l, one has 



= = (-l)'+"n-™(«), 



(45) 



where the overline stands for complex conjugate. The summation in ( HTl) and 

■im=—l' 



below is understood as ^;^o 



Let 



It is well-known that 



h{y) = J2hi{r)Yi{y' 



(46) 



1=0 



(^,m^Z',m')i2(52) — 6ll'6mm', 



(47) 
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where Su' is the Kronecker delta. 
Let 



(48) 



1=0 



where fi := {f,Yi) ^2(^32^ are the Fourier coefficients of /. For sufficiently large 
L one has 

\\f ~ MlHs^) < e/2. (49) 

Thus, if 

1 



/l(/3) + ^ J^e-''^-'h{y)dy 



<e/2, 



(50) 



L2(52) 



then inequality fj40|) holds. Therefore, practically it is sufficient to find /i 
satisfying inequality (ISDjl . Substitute (H^ and ( HTj) into the equation 



e-''^-yh{y)dy = -A7rhi(3), 



B 



and use fllTl) and fH51) to get 



47r(-i)' / r'^ji{kr)hi{r)dr = -47r//^\ < / < L. 



Equation f l52|) can be written as 



(-z)'+M r']i{kr)hi{r)dr = fl"-', V/ > 0. 



Recall that //^^ = //^^^ and hi{r) = /i;,m(r). 
Equation (!53|) has many solutions. 
Denote by h^{r) any function such that 



r"^ ji{kr)h^ {r)dr = 0. 

Then the general solution to equations (1531) has the form 

hi{r) = {-ry^-^fl'^^^.j.^kr) + cihi{r), 0<1<L, 
where q are arbitrarily constants, and 

1 



drr'^jf{kr) 
We have proved the following result. 



(51) 



(52) 



(53) 



(54) 



(55) 



(56) 
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Theorem 3.1 The function ( H6ll . with hi defined in (l55ll /or < I < L and 

hi = for I > L, solves equation (I5T1) . If L = L{e) is sufficiently large, so that 
f H9|) holds, then the function 

L 

h^''\y) = Y,hT\r)Yi{y'). r = \yl y' = ^, (57) 

wt/i h\^\r) defined in fl55l) . satisfies inequality fHOl) . 

Let us give a formula for a potential g such that h = qu approximates /i^^^ 
with a desired accuracy. This g is a solution to the inverse scattering problem 
(IP). 

Given h^^\y) defined in ([57]), let us denote it h{;y) for simplicity. Using 
this function, calculate q{x) formula (fT6l) . If this q G L^{D), then the inverse 
problem (IP) is solved. There is no guarantee that Img < 0. 

If formula f|T6|) does not yield an L'^{D) function, then one uses hs{x) in 
place of h{x) and, as was proved in Theorem 12.41 obtains a potential q5{x) 
which is a solution to (IP). 

Other computational methods can be used for finding h{y) given /(/3). For 
example, one can choose a basis {^j} in L'^{B), B = D is a. ball, and look for 

n 

hn{x) = Y,cfv,{x), (58) 
i=i 

where c^"^ are constants to be found from the minimization problem 

n 

\\f{P)-Y,cf9m\LHs-)=mtn. (59) 
i=i 

Here 

9AP)--=-^J^e-^'^-'vAy)dy- m 

A necessary condition for the minimum in (l59l) is a linear algebraic system for 
the coefficients cj"^ A < j < n. 



13 



References 

[1] S. Agmon, Spectral properties of Schrodinger operators and scattering 
theory, Ann. Scuola Norm. Super. Pisa, 4, (1975), 151-218. 

[2] A.G. Ramm, Recovery of the potential from fixed energy scattering data. 
Inverse Problems, 4, (1988), 877-886. 

[3] A.G. Ramm, Stability of solutions to inverse scattering problems with 
fixed-energy data, Milan Journ of Math., 70, (2002), 97-161. 

[4] A.G. Ramm, Inverse problems, Springer, New York, 2005. 

[5] A.G. Ramm, Uniqueness theorem for inverse scattering problem with non- 
overdetermined data, J.Phys. A, FTG, 43, (2010), 112001. 

[6] A.G. Ramm, Uniqueness of the solution to inverse scattering problem with 
backscattering data, Eurasian Math. Journ (EMJ), 1, N3, (2010), 97-111. 

[7] A.G. Ramm, Distribution of particles which produces a "smart" material. 
Jour. Stat. Phys., 127, N5, (2007), 915-934. 

[8] A.G. Ramm, Wave scattering by many small bodies and creating materials 
with a desired refraction coefficient, Afrika Matematika, 22, Nl, (2011), 
33-55. 

[9] A.G. Ramm, Uniqueness of the solution to inverse scattering problem 
with scattering data at a fixed direction of the incident wave, J. Math. 
Phys., 52, (2011), 123506. 

[10] A.G. Ramm, Scattering by obstacles, D.Reidel, Dordrecht, 1986. 



